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"METODO DI ANALISI DI IMMAGINI RILEVATE DA UN MICRO- 
ARRAY E RELATIVO DISPOSITIVO" 

Campo dell'Invenzione 

La presente invenzione conceme i metodi di analisi automatica di 
immagini e piu in particolare un metodo e relative dispositive di analisi di 
immagini acquisite da micro-array che utilizza una tecnica di filtraggio 
morfologico. 

Background dell'Invenzione 

Le tecniche di analisi di immagini rilevate trovano applicazione in vari 
campi, quali ad esempio Tindividuazione di oggetti fotografati da un satellite 
e Testrazione dei relativi pixel dal resto deH'immagine, 

L' analisi di immagini e di cruciale importanza anche nello studio del 
DNA. In passato e stata sviluppata una tecnica per T analisi del DNA basata 
sui cosiddetti DNA chip. Secondo un metodo sviluppato dalla Stanford 
University, i DNA chip vengono realizzati depositando frammenti di DNA, 
detti anche sonde ("probes"), mediante microdeposizione robotizzata e 
immobilizzati chimicamente su slides da microscopio, ad esempio di 
dimensioni pari a 25mm x 76mm ± 0,5mm e di spessore pari a 1,02mm ± 
0,5mm, di silicio, plastica, nylon o vetro. Le sonde, che corrispondono al 
gene di cui si vuole rilevare la presenza, sono poste sulla slide da 
microscopio secondo una disposizione a matrice (micro-array), in aree 
definite (spot). 

Per verificare se il gene di un individuo contiene una mutazione, 
vengono prelevati dal sangue dell'individuo stesso un campione sospetto di 
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DNA ed un campione che non contiene mutazioni dello stesso gene 
(reference). I campioni di DNA vengono quindi separati in singoli trefoli, 
tagliati in frammenti piu piccoli e maneggevoli ed etichettati con tinte 
fluorescenti. Ad esempio il DNA sano viene colorato con una tinta rossa e 
quello da analizzare con una tinta verde. A questo punto, entrambi i campioni 
sono inseriti nel DNA-chip e, sfinttando lo spontaneo accoppiamento dei 
filamenti complementari (fenomeno delFibridazione), sono lasciati ibridare 
con le sonde contenute nel DNA-chip. 

Se il gene da analizzare non contiene mutazioni, entrambi i campioni 
ibridano alio stesso modo con la sequenza complementare immobilizzata sul 
DNA chip. Se il gene invece contiene mutazioni, il DNA normale (rosso) 
ibrida normalmente con il DNA del chip, ma il DNA da analizzare (verde) 
non presenta ibridazione nella regione corrispondente alia mutazione. 

Utilizzando uno scanner confocale il DNA chip viene quindi sottoposto 
a due diverse scansioni con lunghezza d'onda relativa ai coloranti impiegati. 
Le due immagini ottenute (CY3 e CY5) vengono elaborate da uno specifico 
programma al computer che e in grado di analizzare in base all'intensita della 
fluorescenza se un gene e espresso o meno. 

L'intensita dei pixel componenti gli spot nelle due immagini e 
proporzionale al numero di molecole coloranti, e quindi al numero di sonde 
ibridizzate con il campione contenuto nella posizione corrispondente. 
Confrontando le immagini ottenute (matching) e possibile identificare i geni 
che hanno subito mutazioni. 

Va precisato che qualsiasi analisi di matching e efFettuata solamente 
dopo che il livello del segnale per ogni singolo spot per entrambi i canali e 
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stato normalizzato rispetto al rispettivo tono (livello di grigio) dei pixel di 
background presente. Inoltre e necessaria un'ulteriore operazione di 
normalizzazione tra le immagini ottenute con i due canali dal momento che la 
rispettiva intensita media di spot corrispondenti cambia a seconda del 
colorante usato. 

Con questa tecnica si puo studiare il modo in cui i vari geni sono 
attivati o disattivati per creare certi tipi di cellule. Ad esempio, se un gene e 
espresso in un certo modo, esso genera un cellula normale, mentre se e 
espresso in un altro modo genera una cellula tumorale. Confrontando le 
diverse espressioni dei geni potrebbe essere possibile in fiituro diagnosticare 
con largo anticipo e prevenire la formazione del cancro. Potrebbe inoltre 
essere possibile correlare le variazioni del DNA di una persona con reazioni a 
particolari sostanze, quali le droghe. 

L'elaborazione delle immagini acquisite da micro-array e resa 
complessa dal fatto che le immagini stesse alFintemo di un stesso ciclo di 
esperimenti su un campione di DNA presentano delle caratteristiche 
differenti. Cio e dovuto in parte al processo con cui esse sono generate ed in 
parte aH'acquisizione. Tali operazioni sono influenzate sia dall'intervento 
non ripetibile dell'operatore che dall'imprecisione degli strumenti. 

La Figura 1 mostra 48 spot luminosi di un'immagine, di buona qualita, 
acquisita da micro-array avendo trattato con un colorante (singolo canale) 
campioni di DNA dopo Tibridazione. Filtrando I'immagine e possibile notare 
alcune caratteristiche tipiche di tutte le immagini di micro-array, indicate 
suirimmagine filtrata di Figura 2. 

Gli spot luminosi a sinistra rappresentano "probes" di DNA, visibili in 
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maniera abbastanza nitida neirimmagine filtrata. Questi spot sono tra di loro 
separati, piccoli, compatti e sostanzialmente circolari, localizzati su uno 
sfondo piu scuro (background). Inoltre esiste un mmore (noise) localizzato 
dipendente dal singolo processo di fabbricazione ed imprevedibile. Questo 
rumore causa una variazione del tono (livello di grigio) del background e 
dello spot iuminoso, che rappresenta il segnale utile. 

In generate, I'analisi di immagini di micro-array viene effettuata 
eseguendo le operazioni di: 

• localizzazione (array localization), che consiste nel determinare la 
posizione e la forma degli spot; 

• estrazione (spot extraction), che consiste nell'isoiare singoli spot luminosi, 

• segmentazione (intra-spot segmentation), che consiste nell'esaminare 
ciascuno spot distinguendo pixel di segnale da pixel di sfon 
(background) o di rumore, 

• misura della qualita (spot quality measurement), che consiste nel ricavare 
parametri caratteristici dello spot ed indici che ne determinano la qualita. 

In generale, si puo affermare che le tecniche attualmente utilizzate non 
sono totalmente automatiche ma prevedono un intervento rilevante 
dell'operatore che, pur permettendo di efFettuare un posizionamento (tuning) 
preciso, nei limiti delle caratteristiche del sistema, introduce un errore umano 
non trascurabile. Inoltre I'intervento delFoperatore, non essendo ripetibile 
con precisione, rende i risultati fmali difficilmente confrontabili tra loro. 

II processo di localizzazione e I'operazione che prevede il maggiore 
intervento da parte dell'operatore, che deve centrare ogni singolo spot 
airintemo di una rispettiva cella di una microgriglia. Questa operazione 
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risulta alquanto onerosa dal momento che tipicamente le immagini prese in 
considerazione contengono 10000 e piu spot distribuite su diverse griglie. 

Alio stato deH'arte questa operazione e sostanzialmente 
semiautomatica. Le tecniche di localizzazione permettono di posizionare 
automaticamente la griglia, ma comunque e sempre necessaria una 
sistemazione finale, da parte dell'operatore, per correggere gli errori 
deiralgoritmo di posizionamento. 

Una volta posizionata la griglia, viene generata una mappa binaria che 
defmisce i confini degli spot luminosi sullo sfondo. Questa mappa viene 
usata per isolare gli spot luminosi, che poi vengono esaminati con una tecnica 
di segmentazione. 

Le tecniche di segmentazione piu note per la discriminazione di pixel di 
segnale da quelli di background all' interne di uno spot luminoso sono di 
seguito riportate: 
L Pure Spatial segmentation, 

2. Pure Intensity-based segmentation, 

3. Mixed Spatial/Intensity segmentation, 

4. Mixed Spatial/Statististics segmentation, 

5. Mixed Spatial, Intensity, Statitistics & Morphology segmentation. 

La tecnica di Pure Spatial segmentation si basa suUa assunzione che 
tutti i pixel interni ad un cerchio (o un'altra forma geometrica) di diametro 
preselezionato dall'operatore sono pixel di segnale, mentre quelli interni ad 
una area limitrofa, di forma e distanza selezionata dall'operatore, sono pixel 
di background. In questo caso la discriminazione dei pixel e fatta solamente 
prendendo in considerazione la loro posizione. 
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La tecnica Pure Intensity-based segmentation prevede di considerare 
solamente i pixel alPintemo di una unica area in cui e localizzato lo spot ed 
in base alia intensita del livello dei pixel interni discrimina la parte relativa al 
segnale da quella dello sfondo. In questo caso la discriminazione dei pixel e 
fatta solamente prendendo in considerazione T intensita dei livelli di grigio. 

Con la terza tecnica di segmentazione (Mixed Spatial/Intensity 
segmentation) T intensita dei pixel e Telemento discriminante tra pixel di 
segnale e pixel di sfondo, ma in due diverse regioni, quella circolare e quella 
circostante. Questa tecnica si basa sia su una caratterizzazione spaziale che 
suir intensita dei valori. 

Secondo la tecnica Mixed Spatial/Statististics segmentation si calcola 
con metodi statistici la soglia (livello di grigio) che discrimina un pixel di 
segnale da un pixel di sfondo. L'intensita dei pixel interni alia zona circolare 
e quindi confrontata con tale soglia. 

La tecnica Mixed Spatial, Intensity, Statitistics & Morphology 
segmentation si basa su una conoscenza statistica a priori ottenuta da 
un'analisi locale degli spot, sulla distribuzione delFintensita luminosa e suUe 
caratteristiche morfologiche degli spot stessi. 

Le caratteristiche principali prese in considerazione come indici di 
qualita e di confronto tra spot sono le mediane dei toni (livelli di grigio) 
rispettivamente dei pixels di segnale e dei pixels di sfondo (background). In 
generale non vengono considerate eventuali caratteristiche morfologiche 
degli spot che possono essere utili nella fase di finale di validazione dei 
risultati. 

II fatto che sia necessario Tintervento delFoperatore per analizzare le 
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immagini di micro-array rende I'analisi non riproducibile con sufficiente 
accuratezza. Di conseguenza, qualsiasi operazione di matching e viziata 
aU'origine da un errore umano che puo portare a conclusioni errate. 
SCOPO E SOMMARIO DELL'iNVENZIONE 



immagini rilevate da micro-array in cui I'isolamento degli spot e efFettuato in 
maniera automatica, riducendo inoltre i tempi di elaborazione. 

Secondo un approccio innovativo, gli spot luminosi vengono filtrati e 
isolati senza far ricorso ali'intervento dell'operatore, utilizzando una tecnica 
di filtraggio morfologico. Gli spot cosi isolati vengono successivamente 
analizzati, secondo un altro aspetto innovativo dell'invenzione, con un 
algoritmo fuzzy. 

Piu precisamente, un oggetto delFinvenzione e un metodo di analisi di 
immagini rilevate da micro-array sotto forma di uno o piu spot luminosi su 
uno sfondo, comprendente le seguenti operazioni: 

- determinare forme e relative posizioni di detti spot suU'area del micro- 
array (array localization) generando una mappa binaria di pixels 
definente confini degli spot luminosi sullo sfondo, 

- isolare ciascuno spot mediante un'operazione di estrazione (spot 
extraction) usando la mappa binaria, 

- esaminare gli spot mediante un'operazione di segmentazione (intra-spot 
segmentation) identificando con un criterio prestabilito pixel 
appartenenti ad uno stesso cluster, 

- per ciascuno spot, determinare relativi parametri caratteristici ed indici 
di qualita. 



E stato trovato ed e I'oggetto dell'invenzione un metodo di analisi di 
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II metodo deH'invenzione non richiede Tintervento deH'operatore per 
isolare i singoli spot in quanto la mappa binaria e generata con una tecnica di 
filtraggio morfologico comprendente le operazioni di: 

- filtrare unMmmagine rilevata da microarray con almeno un filtro 
morfologico generando una corrispondente immagine 'marker' dello 
sfondo; 

- ricostruire lo sfondo efFettuando un'operazione di ricostruzione 
suir immagine * marker', generando una corrispondente immagine 
ricostruita di sfondo; 

- generare un'immagine filtrata della luminosita dello sfondo efFettuando 
un'operazione top-hat suH'immagine ricostruita e detta immagine 
rilevata da microarray; 

- generare la mappa binaria mediante un'operazione di thresholding 
suU'immagine filtrata della luminosita dello sfondo. 
Questo metodo e implementato in hardware da un dispositivo di analisi 

di immagini rilevate da micro-array configurato secondo un'architettura di 
rete neurale cellulare. 

Un ulteriore oggetto dell'invenzione e un metodo di identificazione di 
pixel di un'immagine appartenenti ad uno stesso oggetto rappresentato su uno 
sfondo, comprendente scandire i pixel del F immagine, calcolare un valo]/^^ 
caratteristico per ciascun pixel con un algoritmo fuzzy avente coiro ^ 
antecedenti 

- il tono di un pixel, 

- la distanza tra il tono del pixel e il valore medio del tono di rispettivi 
pixel di sfondo, 
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- il quadrato di questa distanza; 

e contemplante le operazioni di calcolare il valore medio del tono di pixel di 
sfondo e giudicare i pixel scanditi come appartenenti ad uno stesso oggetto se 
il loro valore caratteristico eccede una soglia prestabilita. 

Questo metodo e implementato in un relativo sistema fuzzy di 
identificazione. 

L'invenzione e piu precisamente defmita nelle annesse rivendicazioni. 
Breve Descrizione dei Disegni 

I diversi aspetti e vantaggi dell'invenzione risulteranno ancor piu 
evidenti attraverso una descrizione dettagliata dell'invenzione facendo 
riferimento ai disegni allegati, in cui: 

la Figura 1 mostra un'immagine rilevata da micro-array; 

la Figura 2 mostra Timmagine della Figura 1 filtrata; 

la Figura 3 mostra un diagramma a blocchi di una forma preferita di 
realizzazione del metodo di analisi di immagini deH'invenzione; 

la Figura 4 mostra lo schema dei filtraggi morfologici da effettuare su 
un'immagine rilevata da micro-array secondo una forma preferita di 
realizzazione del metodo di analisi di immagini dell'invenzione; 

la Figura 5 mostra uno spot luminoso da filtrare, con la relativa mappa 
dei livelli; 

le Figure da 6 a 9 mostrano spot, con le relative mappe dei livelli, 
ottenuti mediante successivi filtraggi dello spot di Figura 5 con i filtri 
morfologici indicati in Figura 4; 

la Figura 10 mostra uno spot, con la relativa mappa dei livelli, di 
un'immagine di sfondo ottenuta efFettuando un'operazione di ricostruzione 
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sullo spot di Figura 9; 

la Figura 11 mostra uno spot filtrato da rumore, con la relativa mappa 
dei livelli, ottenuto efFettuando un'operazione di top-hat sugli spot di Figure 
5 e 10; 

la Figura 12 e un diagramma di flusso deU'operazione di array 
localization secondo una forma di realizzazione preferita del metodo di 
analisi di immagini deirinvenzione; 

la Figura 13 descrive schematicamente I'operazione di estrazione di 
spot del metodo di analisi di immagini dell'invenzione; 

la Figura 14 e un diagramma di flusso di un algoritmo di 
identificazione di cluster (clustering) di un'immagine; 

le Figure da 15 a 20 illustrano il funzionamento delFalgoritmo di 
identificazione di cluster (clustering) descritto in Figura 14; 

le Figure da 21 a 25 illustrano Talgoritmo di identificazione di cluster 
(clustering) descritto in Figura 14 per eseguire Toperazione di estrazione di 
spot di un'immagine rilevata da micro-array; 

la Figura 26 descrive un sistema di estrazione (spot extraction) di spot 
luminosi su uno sfondo di immagini rilevate da micro-array; 

la Figura 27 e un diagramma di simulazione del sistema fuzzy 
deirinvenzione in corrispondenza di pattern di allenamento corrispondenti a 
due spot con livelli differenti; 

la Figura 28 e un diagramma di fiinzioni membro usate nell' algoritmo 
fuzzy deirinvenzione relative al tono di un pixel; 

la Figura 29 e un diagramma di funzioni membro usate nell' algoritmo 
fiizzy deirinvenzione relative alia distanza tra il tono di un pixel e il valore 
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medio del tono dei rispettivi pixel di sfondo; 

la Figura 30 e un diagramma di fiinzioni membro usate neiralgoritmo 
fuzzy dell'invenzione relative al quadrato della distanza tra il tono di un pixel 
e il valore medio del tono dei rispettivi pixel di sfondo; 

la Figura 31 e un diagramma dei conseguenti deiralgoritmo fixzzy 
deirinvenzione; 

le Figure da 32 a 40 mostrano come calcolare parametri caratteristici 
ed indici di qualita degli spot isolati con il metodo deirinvenzione. 

Descrizione di Alcune Forme di Realizzazione dell'Invenzione 
Secondo un aspetto innovativo deirinvenzione, Toperazione di 
localizzazione e efFettuata mediante filtraggi morfologici. Poiche le forme di 
interesse degli spot sono approssimabili a dei cerchi di diametro non 
superiore ad una certa lunghezza, e quindi appartenenti ad una particolare 
classe di figure geometriche simmetriche rispetto al lore centre, e possibile 
progettare dei filtri morfologici servendosi di opportune Opening direzionali. 

Preferibilmente, si usa una cascata di quattro filtri morfologici con 
rispettive Opening che impiegano come elemento strutturante (structuring 
set) segmenti di lunghezza non superiore al massimo diametro dello spot e 
diretti lungo le due direzioni fondamentali di un piano cartesiano (X, Y) e 
lungo le due diagonali principali (Di, D2), come schematicamente riassunto 
in Figura 4. 

In alternativa, e possibile utilizzare anche filtri morfologici diversi da 
quelli citati come ad esempio un unico filtro morfologico avente uno 
structuring set formato da quattro segmenti che si intersecano in un unico 
punto e suddividono il piano in otto angoli uguali. Da prove effettuate si e 
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visto che in quest'ultimo caso si ottengono, in genere, risultati peggiori 
rispetto al filtraggio eseguito con i quattro filtri morfologici citati. 

Filtrando lo spot di Figura 5 in successione con i filtri morfologici di 
Figura 4, si ottengono rispettivamente gli spot delle Figure da 6 a 9. Gli spot 
filtrati con i quattro filtri morfologici, come quello rappresentato in Figura 9, 
costituiscono la cosiddetta immagine 'marker' dello sfondo. 

Successivamente si effettua un'operazione di ricostruzione 
(Reconstmctiori) sulPimmagine 'marker' ottenendo un'immagine ricostruita 
di sfondo costituita da spot come, quello di Figura 10. Infine si esegue 
un'operazione top-hat con Timmagine originaria e Timmagine ricostruita di 
sfondo, ottenendo un'immagine filtrata della luminosita dello sfondo 
costituita da spot come quello di Figura 1 1, corrispondente alio spot di Figura 
5, in cui Tintensita dello sfondo e stata praticamente annullata. 

Infine viene generata la rispettiva mappa binaria mediante 
un'operazione di thresholding, che sostanzialmente consiste nel confrontare 
con una soglia prestabilita la mappa dei livelli dell' immagine ottenuta con 
I'operazione di top-hat. se il livello di grigio di un pixel eccede la soglia 
prestabilita, viene generato un corrispondente pixel della mappa binaria 
avente un valore logico attivo (1), altrimenti viene generato un pixel della 
mappa binaria con valore logico nullo (0). 

A causa di effetti di bordo, ma anche poiche spesso vengono elimi^S| 
completamente pixel connessi appartenenti al background o ad elemenfi^d?^^ 
rumore, I'operazione di reconstruction non permette di ottener^^, 
perfettamente 1' immagine dello sfondo degli spot. L'operazione top-hat 
fomisce quindi un'immagine in toni di grigio in cui oltre ad avere gli spot e 
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presente anche un rumore impulsivo. Di conseguenza, anche la mappa binaria 
sara corrotta da rumore. 

La distribuzione del rumore impulsivo che corrompe la mappa binaria e 
tale da consentirne una facile rimozione impiegando due operazioni di 
erosion, con elemento strutturante costituito da due maschere circolari di 
raggio diverse. 

Per motivi che saranno piu evidenti in seguito, risulta desiderabile 
ottenere una mappa binaria che consenta di delimitare insieme alio spot 
luminoso, anche un suo intomo di ampiezza sufficiente a calcolare con una 
certa afFidabilita il valore locale di sfondo. A tale scopo, secondo una forma 
di realizzazione del metodo dell'invenzione si effettua un'operazione di 
dilation mediante una maschera circolare di diametro superiore al diametro 
dello spot. Questa operazione di dilation consente di compensare gli efFetti 
indesiderati delle operazioni di erosion, che sostanzialmente consistono nel 
restringimento (sulla mappa binaria) della zona dello spot luminoso che non 
contiene solo pixel di sfondo o di rumore. 

Un diagramma di flusso che riassume gli step di processo di 
un'operazione di localizzazione (array localization) secondo il metodo di 
analisi dell'invenzione, e mostrato inFigura 12. 

Convenientemente, I'operazione di localizzazione (array localization) 
sara eseguita da un sistema di array localization organizzato secondo 
un'architettura di rete neurale cellulare del tipo divulgate in [1]. 

L'immagine ottenuta con il processo di array localization e una mappa 
logica binaria, esemplificata in Figura 13, contenente la sagoma di ogni spot. 



[I] EP 1 182 602 di Arena P. et al., a nome deUa stessa richiedente. 
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I pixel in nero hanno valore logico attivo (1) ed indicano la presenza di uno 
spot, mentre i pixel in bianco hanno valore logico nuUo (0) e corrispondono a 
pixel di sfondo o di rumore. 

Al fine di estrarre gli spot e possibile usare un qualsiasi algoritmo noto 
di identificazione di cluster (clustering) da immagini basato su proprieta dei 
pixels. 

Un algoritmo di identificazione di cluster particolarmente adatto alio 
scopo e Talgoritmo divulgato in [2], implementato utilizzando la distanza 
euclidea come metrica e un criterio a soglia (clustering condition) per 
individuare pixel di uno stesso cluster. 

Questo algoritmo, che verra descritto in seguito, puo essere applicato 
anche per eseguire operazioni diverse da quella di estrazione di spot. Esso e 
molto conveniente perche flinziona qualunque sia il predicato specificate 
dali'utente per distinguere pixel di segnale da pixel di sfondo, senza che la 
successione degli step di processo sia modificata. Inoltre esso permette di 
utilizzare predicati basati su "soft computing" e ha il vantaggio di non 
richiedere rintervento delFutente. Infine esso non richiede alcuna condizione 
iniziale (seme dell' algoritmo), a difFerenza di tanti altri algoritmi che 
producono risultati diversi a seconda del seme inizialmente scelto. 

Un diagramma di flusso generate delF algoritmo di identificazione 
citato, utilizzato per eseguire Toperazione di estrazione del metodo di analisi 
di immagini dell'invenzione, e rappresentato in Figura 14, 

II primo passo GENERAL Clustering consiste nello scandire per 
colonna (o per riga) i pixel deH'immagine filtrata, raggruppando in cluster 

[2] AJessi Enrico, Coco Salvatore, Pappalardo Giuseppe, Capizzi Giacomo, "A NEW 

CLUSTERING BASED SYSTEM FOR AUTOMATED OBJECT RECOGNITION'' - 
SOCO '98. - Universita di Catania. 
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elementari i pixel di una stessa colonna (riga) scanditi consecutivamente, cioe 
pixel adiacenti della stessa colonna (riga), che soddisfano un criterio 
prestabilito. Ciascun cluster elementare e definito da un numero identificativo 
INDEX, da un numero di colonna y (riga) e da coordinate minima Xmin e 
massima Xmax del cluster sulla colonna (riga) 

II criterio particolare usato per eseguire I'operazione di estrazione di 
spot (spot extraction), consiste nel raggruppare in uno stesso cluster 
elementare pixel scanditi consecutivamente che appartengano alia stessa 
colonna (riga) e che corrispondano a pixel della mappa binaria aventi lo 
stesso valore logico. 

Per meglio comprendere il funzionamento delFalgoritmo, se ne dara 
una breve descrizione generale facendo riferimento alle Figure da 15 a 20. 

LaFigura 15 mostra colonne di pixel di un'immagine gia raggruppati in 
cluster elementari in seguito all'operazione di GENERAL CLUSTERING 
eseguita con un certo criterio. Ad esempio, questo criterio potrebbe essere 
determinato da un predicato basato su "soft computing", oppure potrebbe 
essere un criterio basato su soglia. II numero identificativo del cluster 
elementari e riportato a sinistra di ciascun cluster. 

I pixel di uno stesso cluster elementare sono rappresentati come pixel di 
un'immagine binaria aventi valore logico attivo (1), mentre i pixel di sfondo 
o di rumore sono rappresentati come pixel aventi valore logico nullo (0), 

I cluster elementari sono processati con una prima operazione di 
fusione MERGE 1, efFettuata in due passi: prima si individua, per ogni cluster 
elementare c in una certa colonna (riga), un insieme di cluster elementari S 
nella colonna (riga) immediatamente precedente che confmano con tale 
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cluster elementare; poi in questo insieme S si identifica un cluster "vincitore" 
w che ha il maggior numero di pixel di confine con il cluster c e si uguaglia il 
sue numero identificativo con quello del cluster c. II risultato dell'operazione 
MergeI suirimmagine binaria di Figura 15 e mostrato in Figura 16. 

I cluster deir insieme S diversi dal cluster vincitore vengono individuati 
come cluster "perdenti" e ad essi e assegnata una rispettiva tripla 
comprendente il loro numero identificativo, il numero identificativo del 
rispettivo cluster vincitore e il numero di cplonna cui appartengono (BOUND 
Reporting, Report Processing). 

I numeri identificativi dei cluster perdenti vengono cambiati nei numeri 
identificativi dei rispettivi cluster vincitori mediante I'operazione MERGE2, il 
cui risultato e mostrato in Figura 17. A questo punto i cluster elementari 
appartenenti ad uno stesso oggetto sono individuati da uno stesso numero 
identificativo. 

Per comodita si possono assegnare ai diversi oggetti rappresentati 
numeri identificativi consecutivi con Toperazione di CLUSTER Sorting, 
come mostrato in Figura 18. 

Infine con I'operazione CLUSTER Extraction, si estraggono 
dairimmagine originaria di Figura 19, associata alPimmagine binaria di 
Figura 15, i pixel corrispondenti a cluster elementari aventi lo stesso numero 
identificativo, isolando cosi i singoli oggetti rappresentati, come mostrato 
Figura 20. 



Applicando Talgoritmo appena descritto ad una mappa binaria ottenuta^ 
mediante I'operazione di array localization e alia relativa immagine filtrata, i 
singoli spot luminosi sono isolati senza I'intervento dell'operatore, come 
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descritto dalle Figure da 21 a 25. In Figura 21 e mostrata in forma stilizzata 
una mappa binaria. La Figura 22 mostra il risultato deH'operazione MERGEl 
(roperazione MERGE2 nel caso di Figura 21 non cambia nulla) e la Figura 23 
il risultato delPoperazione di CLUSTER SORTING. Le Figure 24 e 25, 
analogamente alle Figure 19 e 20, illustrano Testrazione degli spot luminosi 
dall'immagine rilevata da microarray. 

Un possibile sistema di estrazione (spot extraction) di spot luminosi su 
uno sfondo di immagini rilevate da micro-array e mostrato in Figura 26. Esso 
comprende un sottosistema di scansione Extractor di pixel di 
un'immagine, un sottosistema di determinazione (FEATURES EXTRACTION, 
Clustering Condition) di cluster elementari e un sottosistema di 
elaborazione CLUSTERING, II blocco circuitale FEATURES EXTRACTION ricava 
dai pixel scanditi il rispettivo livello di grigio e fornisce questa informazione 
al blocco in cascata CLUSTERING Condition, che raggruppa i pixel in cluster 
elementari. Infme il sottosistema di elaborazione CLUSTERING dei cluster 
elementari produce in uscita pixels di cluster CLUSTER 1, CLUSTER N 
rappresentanti gli spot luminosi secondo Toperazione di estrazione prima 
descritta. 

Gli spot luminosi cosi isolati dalFimmagine vengono processati 
singolarmente con I'operazione di segmentazione (intra-spot segmentation), 
al fine di individuare oggetti in essi rappresentati, 

L'operazione di segmentazione sostanzialmente consiste 
nel r individuare pixel appartenenti ad uno stesso oggetto usando un criterio di 
similitudine tra pixel. Teoricamente e possibile utilizzare difFerenti criteri di 
similitudine basati suU'intensita del pixel e su diverse proprieta deU'intomo 
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stesso del pixel. Owiamente la scelta del criterio deve essere efFettuata 
considerando il problema specifico da risolvere. 

Se, per esempio, si volessero estrarre tutti gli oggetti di una immagine 
costituiti dai contorai degli oggetti rappresentati nell' immagine stessa e se 
rimmagine fosse priva di rumore, si potrebbe basare il criterio solamente sul 
valore del gradiente di una fmestra di dimensione 3x3 o 5x5 attomo al pixel 
considerate. Nel caso in cui 1' immagine e corrotta da rumore, non basta 
applicare un predicato semplice come quello appena menzionato ma bisogna 
utilizzare dei criteri basati su diverse proprieta, molto piu onerosi dal punto di 
vista computazionale. 

II ragionamento umano per I'individuazione di un oggetto aU'interno di 
un' immagine e molto complesso, basato su molteplici parametri e a volte 
difficilmente definibile dal punto di vista matematico. Proprio per questo 
motivo, per risolvere i problemi piu complessi di identificazione vengono 
utilizzati criteri defmiti da predicati basati su soft computing, in alternativa ai 
criteri definiti da predicati a soglia. 

In particolare si pud usare un predicato neurale di tipo MLP (Multi 
layer perceptron) ed un altro basato suUa logica fiizzy. Entrambi i predicati 
permettono di superare i problemi che si pongono in situazioni che 
richiederebbero criteri difficilmente defmibili matematicamente; inoltre 
permettono una decisione multicriterio e conferiscono al sistema stesso 
robustezza qualora si vogliano individuafe oggetti in immagini corrotte da 
rumore relativamente elevato. 

Secondo un altro aspetto innovativo della presente invenzione, 
I'operazione di segmentazione (intra-spot segmentation) viene efFettuata 
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mediante un algoritmo fuzzy, che puo essere definito da regole linguistiche 
determinate direttamente da un esperto oppure definito a partire da un file di 
learning che contiene i pattern delle proprieta e le corrispondenti uscite in 
relazione al problema da risolvere. 

La sagoma dello spot della mappa binaria ottenuta con Toperazione di 
array localization e piu grande di quella dello spot reale, in modo da 
selezionare con lo spot anche la zona circostante di sfondo (background). II 
tono (livello di grigio) dello sfondo nelle zone limitrofe dello spot isolato 
serve per il calcolo di un fattore di normalizzazione e pertanto deve essere 
determinato con precisione. 

Questo fattore di normalizzazione e importante perche i livelli di grigio 
di immagini ottenute quando si analizzano campioni di DNA dipendono 
anche dal tipo di colorante usato per evidenziare I'ibridazione. Per poter 
confrontare (matching) quindi immagini ottenute con coloranti diversi, i toni 
degli spot luminosi devono essere prima normalizzati rispetto al tono dello 
sfondo. Con altrettanta cura deve essere determinato il tono di ciascun pixel 
soprattutto in presenza di spot facilmente confondibili con lo sfondo 
(background) locale. 

E importante ricordare che pixel di background possono trovarsi 
annidati anche all'interno dello spot. Questo fa si che i criteri da utilizzare per 
distinguere pixel di segnale da pixel di sfondo all'interno di uno spot in 
genere sono complicati o diflficilmente formalizzabili. Per owiare a questa 
difficolta e stato ideato un criterio fiizzy per distinguere pixel di segnale da 
pixel di sfondo o di rumore airintemo di uno spot luminoso. Secondo tale 
criterio fuzzy, si scandiscono pixel di uno spot e si calcola un valore 
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caratteristico per ogni pixel scandito: se tale valore caratteristico eccede una 
soglia prestabilita, allora il pixel viene giudicato come appartenente ad un 
oggetto (cluster) rappresentato nello spot. 

La scelta degli antecedenti impiegati neiralgoritmo fiizzy e stata 
efFettuata dopo aver studiato ristogramma dei livelli di grigio degli spot. 

E stato notato che ogni spot ha un istogramma normalmente trimodale 
o bimodale a seconda che esso e o non e corrotto da rumore. In particolare 
mentre il primo e il secondo modo corrispondono rispettivamente al 
background e al segnale, il terzo modo e relativo al rumore. Sono cosi stati 
scelti i seguenti 3 antecedenti per un algoritmo fuzzy : 

1 . il tono del singolo pixel; 

2. la distanza (5) del livello di grigio del pixel considerato dal valor medio 
del tono di pixel sfondo; 

3. il quadrato di questa distanza (5^). 

Le regole, la forma ed il posizionamento delle fiinzioni membro 
(membership functions) sono state generate automaticamente servendosi di 
una rete neurale. 

Preferibilmente, ma non obbligatoriamente, airinterno dello spot sono 
identificate preliminarmente due zone: la zona True signal contenente pixel 
di segnale e la zonei False signal contenente pixel di rumore o di sfondo. Per 
la separazione dei pixel dello spot in una zona contenente di pixel di se^ra 



ed in una zona contenente di pixel di sfondo e utilizzato Talgoritmo wz^y 

W 

con i 3 antecedenti sopra citati, ciascun antecedente avente 3 funzicx 
membro (membership functions). 

Per I'apprendimento (learning) di un sistema fuzzy destinato ad 
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implementare ralgoritmo, sono stati individuati gli spot piu rappresentativi 
della maggior parte dei casi che si possono verificare. Ad ogni pattern e 
associate un valore logico 1 o 0 a seconda che il pixel considerato sia un 
pixel di segnale o di sfondo. 

In Figura 27 e mostrato un diagramma di simulazione deiralgoritmo 
fuzzy in corrispondenza dei pattern di allenamento corrispondenti a due spot 
con livelli di grigio completamente difFerenti. Le Figure da 28 a 30 mostrano 
le fiinzioni membro delle variabili di ingresso delPalgoritmo fuzzy 
rispettivamente del tono di grigio, della distanza 5 e della distanza al quadrato 
6^. I conseguenti delPalgoritmo fiizzy sono riportati in Figura 31. 

Maggiori dettagli sul predicate fuzzy deiralgoritmo di segmentazione 
(intra-spot segmentation) sono esposti nel seguente rapporto: 



SISTEMA FUZZY CON FUNZIONI MEMBRO GAUSSIANE 

\x-m) 



Formula della fiinzione membro (shape): y ~ exp 
in cui: 

X e il valore d' ingresso; 
/w e la media; 
5 e la varianza. 
Nel seguito sara indicata con GAUSS(x, m, s). 



*** Variabili d'ingresso *** 
Variabile d' ingresso 1: tono 

Fuzzy Set Name: Centroid: 

SET1_1 25.182846 

Shape: GAUSS(xi, 25.182846, 9.731716) 

Fuzzy Set Name: Centroid: 

SET1_2 22.915066 

Shape: GAUSS(xi, 22.915066, 1.789337) 

Fuzzy Set Name: Centroid: 

SET1_3 62.334343 

Shape: GAUSS(xi, 62.334343, 30.384819) 



2s' 



J 



Variabile d'ingresso 2: distanza 
Fuzzy Set Name: 
SET2 1 1.411060 



Centroid: 
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Shape: GAUSS(x2, 1.411060, 1.927289) 

Fuzzy Set Name: Centroid: 

SET2_2 11.101624 

Shape: GAUSS(x2, .l 1.101624, 7.363041) 

Fuzzy Set Name: Centroid: 

SET2_3 7.316561 

Shape: GAUSS(x2, 7.316561, 31.509554) 

Variabile d'ingresso 3: distanza al quadrato 
Fuzzy Set Name: Centroid: 
SET3_1 88.555168 

Shape: GAUSS(X3, 88.555168, 1090.161377) 

Fuzzy Set Name: Centroid: 

SET3_2 1176.142700 

Shape: GAUSS(x3, 1176.142700, 310.110229) 

Fuzzy Set Name: Centroid: 

SET3_3 1488.746460 

Shape: GAUSS(x3, 1488.746460, 308.677734) 

*** Variabili di uscita *** 

Variabile di uscita 1: CONSl 

Fuzzy Set Names: Centroids: 

CONS1_l 0.006696 

CONSl_17 -0.413813 

CONSl_21 1.986579 

CONSl_63 1.001855 

CONSl 64 0.999214 



II criterio fuzzy appena descritto puo convenientemente essere 
utiiizzato per eseguire lo step GENERAL CLUSTERING deiralgoritmo di 
indentificazione divulgate in [2] per eseguire I'operazione di segmentazione. 
Secondo un'implementazione preferita del metodo di analisi di immagini 
deirinvenzione, Toperazione di segmentazione e eseguita efFettuando i passi 
di processo schematizzati in Figura 14 utilizzando pero questo criterio fiizzy 
invece del criterio basato su soglia, che e stato usato per estrarre i singoli 
spot. 

Owiamente, un sistema per eseguire I'operazione di segmentazione 
(intra-spot segmentation) sara come quello riportato in Figura 26 e descritto 
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con riferimento all'operazione di estrazione di spot. In questo caso perd il 
blocco Features Extraction ricava per ogni pixel scandito i rispettivi tre 
antecedent! del criterio fiizzy (tono, 6, 5^), e fomisce questa informazione al 
blocco Clustering Condition, il quale raggruppa i pixel in cluster 
elementari usando il criterio fuzzy. II blocco CLUSTERING esegue nel modo 
gia descritto gli altri passi deiralgoritmo. 

I cluster di uscita CLUSTER 1, CLUSTER N rappresentano in questo 
caso oggetti (costituiti da pixel di segnale) rappresentati nello spot esaminato. 

Una volta completata I'operazione di segmentazione (intra-spot 
segmentation), ogni singolo spot del micro-array viene infine esaminato per 
ottenere dei parametri caratteristici e indici di qualita. Questi parametri 
permettono di caratterizzare uno spot ed eventualmente di efFettuare anche 
un'analisi comparativa tra le immagini ottenute usando i coloranti Cy3 e Cy5 
per evidenziare Tibridazione in campioni di DNA. 

La seguente tabella illustra alcuni parametri utili per caratterizzare uno 
spot e il modo in cui essi sono calcolati, facendo riferimento alia figura citata. 



Parametro 


Significato 


Centro dello spot 


Coordinate del baricentro dello spot 
Figura 32 


Media del segnale 


Media dei livelli di grigio dei pixel dello spot 
appartenenti alia classe True signal 

Figura 33 


Media del background 


Media dei livelli di grigio dei pixel dello spot 
appartenenti al perimetro della shape 
Figura 34 


Mediana del segnale 


Mediana dei livelli di grigio dei pixel dello spot 
appartenenti classe True signal 
Figura 35 


Mediana del background 


Mediana dei livelli di grigio dei pixel dello spot 
appartenenti al perimetro della shape 
Figura 36 


FOMV 


Media dell' uscita del sistema fuzzy per i pixel 
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assegnati alia classe True signal. 


Fattore di simmetria 


Rapporto tra altezza e larghezza del minimo 
rettangolo che include lo spot 
Figura 37 


Area dello spot 


Numero di pixel dello spot 
Figura 38 


Area del background 


Numero di pixel del perimetro della shape 
Figura 39 


Area del segnale 


Numero di pixel appartenenti alia classe True 

signal 
Figura 40 


Fattore di correzione 


DifFerenza tra la mediana del livello di grigio dei 
pixel del segnale e la mediana del livello di grigio 
dei pixel del background 



Un parametro importante per la valutazione della qualita dello spot 
estratto e quelle indicate cen la sigla FOMV, date dalla media dell'uscita 
fiizzy per pixel di segnale, che e un indice della qualita della selezione dei 
pixel di segnale dai pixel di sfondo. 

Quanto piu I'uscita fuzzy si approssima ad 1, tante piii si puo 
considerare affidabile I'individuazione di un pixel come pixel di segnale e 
non di sfondo. Facendo una media dell'uscita fuzzy, si ottiene quindi un 
indice della qualita della selezione dei pixel di segnale tra i pixel di sfondo. 

METODO di IDENTIFICAZIONE di pixel appartenenti ad UNO STESSO 
CLUSTER 

Un ulteriere oggetto della presente invenzione e un metodo di 
identificazione di pixel di un'immagine appartenenti ad uno stesso cluster. 

Tutti i metodi di identificazione di pixel di uno stesso oggetto 
rappresentato si basano su un criterio di similitudine tra pixel. E possibile 
utilizzare differenti criteri di similitudine basati sulPintensita del pixel e su 
diverse proprieta deirintomo stesso del pixel. Owiamente la scelta del 
criterio deve essere efFettuata considerando il problema specifico da 
risolvere. 
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Se, per esempio, si vogliono identificare gli oggetti di una immagine 
aventi contomi definiti e se 1' immagine fosse priva di rumore, il criterio di 
similitudine potrebbe essere basato solamente sul valore del gradiente di una 
finestra di dimensione 3x3 o 5x5 attomo al pixel considerate. Nel caso in cui 
rimmagine e corrotta da rumore, non basta applicare un predicate semplice 
come queilo appena menzionato ma bisogna utilizzare altri criteri molto piu 
onerosi dal punto di vista computazionale. 

II ragionamento umano per Tindividuazione di un oggetto airinterno di 
un' immagine e molto complesso, basato su molteplici parametri e a volte e 
difficilmente defmibile dal punto di vista matematico. Proprio per questo 
motivo, per risolvere i problemi piu complessi di identificazione, vengono 
utilizzati criteri definiti da predicati basati su soft computing, in alternativa ai 
criteri definiti da predicati a soglia. 

In particolare si puo usare un predicate neurale di tipo MLP (Multi 
layer perceptron) ed un altro basato su un algoritmo fuzzy, che puo essere 
definite da regole linguistiche determinate direttamente da un esperte oppure 
definite a partire da un file di learning che contiene i pattern delle proprieta e 
le corrispondenti uscite in relazione al problema da risolvere. Entrambi i 
predicati permettono di superare i problemi che si pengone in situcizioni che 
richiederebbero criteri difficilmente definibili matematicamente; inoltre 
permettono una decisione multicriterie e conferiscono al sistema stesso 
robustezza qualora si vegliane individuare oggetti in immagini corrette da 
rumore relativamente elevate. 

II metedo di identificazione delFinvenzione si basa su un criterio fuzzy 
per distinguere pixel di un oggetto da identificare da pixel di sfondo o di 
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aimore. Uno scopo deirinvenzione e quelle di fomire un metodo di 
identificazione di pixel appartenenti ad uno stesso oggetto che possa essere 
usato per eseguire i'operazione di segmentazione (intra-spot segmentation) di 
spot luminosi su uno sfondo all' interne di un metodo di analisi di immagini 
rilevate da micro-array. 

Secondo tale metodo dell'invenzione, si scandiscono pixel di immagine 
e si calcola con un algoritmo fuzzy un valore caratteristico per ogni pixel 
scandito: se tale valore caratteristico eccede una soglia prestabilita, allora il 
pixel viene giudicato come appartenente ad un oggetto (cluster) 
rappresentato. 

La scelta degli antecedenti impiegati nell' algoritmo fuzzy e stata 
effettuata dopo aver studiato I'istogramma dei livelli di grigio di tipici spot 
luminosi di immagini che si ottengono nello studio del DNA. 

E state notato che ogni spot ha un istogramma normalmente trimodale 
o bimodale a seconda che esso e o nen e corrotto da rumore. In particolare 
mentre il prime e il secondo mode corrispondene rispettivamente al 
background e al segnale, il terzo modo e relative al rumore. Sone cosi stati 
scelti i seguenti 3 antecedenti per un algoritmo fuzzy: 

1 . il tone del singolo pixel; 

2. la distanza (6) del livello di grigio del pixel considerate dal valor medio 
del tone di pixel sfondo; 

3. il quadrate di questa distanza (5^). 

Sostanzialmente, il metodo dell'invenzione consiste nello scandire 
pixel deir immagine, calcolare un valore caratteristico per ciascun pixel con 
un algoritmo fiizzy i cui antecedenti sone 
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• il tono di un pixel, 

• la distanza tra detto tono e il valore medio del tono dei rispettivi pixel 
di sfondo, 

• il quadrato di detta distanza; 

e contemplante le operazioni di calcolare il valore medio del tono di pixel di 
sfondo, giudicando i pixel scanditi come appartenenti ad uno stesso oggetto 
se il rispettivo valore caratteristico eccede una soglia prestabilita. 

Preferibilmente, ciascun antecedente ha tre funzioni membro 
(membership functions), ciascuna di esse essendo una fiinzione gaussiana con 
media e varianza predeterminate. 

Le Figure da 28 a 30 mostrano alcune funzioni membro delle variabili 
di ingresso deH'algoritmo fuzzy rispettivamente del tono di grigio, della 
distanza 5 e della distanza al quadrato 5^. I conseguenti deiralgoritmo fiizzy 
sono riportati in Figura 31. Maggiori dettagli sul criterio fuzzy descritto 
facendo riferimento alle Figure da 28 a 31 sono esposti nel seguente rapporto: 



2s' 



J 



SISTEMA FUZZY CON FUNZIONI MEMBRO GAUSSIANE 
Formula della fiinzione membro (shape): y = exp 

V 

in cui: 

jc e il valore d' ingresso; 
777 e la media; 
5 e la varianza. 
Nel seguito sara indicata con GAUSS(x, w, 5). 

*** Variabili d' ingresso *** 

Variabile d' ingresso 1 : tono 
Fuzzy Set Name: Centroid: 
SET1_1 25.182846 
Shape: GAUSS(jci, 25,182846, 9,731716) 
Fuzzy Set Name: Centroid: 
SET1_2 22.915066 
Shape: GAUSS(xi, 22,915066, 1.789337) 
Fuzzy Set Name: Centroid: 
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SET1_3 62.334343 

Shape: GAUSS(xi, 62.334343, 30.384819) 

Variabile d'ingresso 2: distanza 

Fuzzy Set Name: Centroid: 

SET2_1 1.411060 

Shape: GAUSS(x2, 1.411060, 1.927289) 

Fuzzy Set Name: Centroid: 

SET2_2 11.101624 

Shape: GAUSS(x2, 11.101624, 7.363041) 

Fuzzy Set Name: Centroid: 

SET2_3 7.316561 

Shape: GAUSS(x2, 7.316561, 31.509554) 

Variabile d'ingresso 3: distanza al quadrato 
Fuzzy Set Name: Centroid: 
SET3_1 88.555168 

Shape: GAUSS(x3, 88.555168, 1090.161377) 

Fuzzy Set Name: Centroid: 

SET3_2 1176.142700 

Shape: GAUSS(X3, 1176.142700, 310.110229) 

Fuzzy Set Name: Centroid: 

SET3_3 1488.746460 

Shape: GAUSS(x3, 1488.746460, 308.677734) 

*** Variabili di uscita *** 



Variabile di uscita I: CONSl 

Fuzzy Set Names: 

CONSl_l 

C0NS1_17 

C0NS1_21 

CONSl_63 

CONSl 64 



Centroids: 

0.006696 

-0.413813 

1.986579 

1.001855 

0.999214 
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Secondo una forma preferita di realizzazione della presente invenzione, 

il criterio fuzzy appena descritto puo essere convenientemente utilizzato 

neiralgoritmo di estrazione divulgate in [2] per estrarre pixel di oggetti 

rappresentati in un'immagine. Un diagramma di flusso di questj^^ 

If 

realizzazione preferita deiralgoritmo di identificazione dell'invenzione & 
rappresentato in Figura 14. 

II primo passo GENERAL Clustering consiste nello scandire per 



ivi^HeA^i>AvB<C>l;;IK>l 



29 



N 



CAMERA Dl COMMERCIO 
INDUSTRIA ARTIGIANATO AGRICOLTURA 
Ufficio Brevetti e Marchi 
VARESE 



2 2 nOV. 2002 



VA/ 2002/A/ 0 0 6 0 




colonna (o per riga) i pixel deirimmagine filtrata, raggaippando in cluster 
elementari i pixel di una stessa colonna (riga) scanditi consecutivamente, cioe 
pixel adiacenti della stessa colonna (riga), che secondo il criterio fiizzy 
descritto sono giudicati come appartenenti ad uno stesso oggetto. Ciascun 
cluster elementare e definito da un numero identificativo INDEX, da un 
numero di colonna y (riga) e da coordinate minima Xmin e massima Xmax del 
cluster sulla colonna (riga)>'. 

Per meglio comprendere questa realizzazione del metodo di 
identificazione deirinvenzione, si faccia riferimento alle Figure da 15 a 20. 

La Figura 15 mostra colonne di pixel di un'immagine gia raggruppati in 
cluster elementari in seguito al passo di processo di GENERAL CLUSTERING. II 
numero identificativo dei cluster elementari e riportato a sinistra di ciascun 
cluster. 

I pixel di uno stesso cluster elementare sono rappresentati come pixel di 
un'immagine binaria aventi valore logico attivo (1), mentre i pixel di sfondo 
o di rumore sono rappresentati come pixel aventi valore logico nullo (0). 

I cluster elementari sono processati con una prima operazione di 
fusione MERGEl, efFettuata in due passi: prima si individua, per ogni cluster 
elementare c in una certa colonna (riga), un insieme di cluster elementari S 
nella colonna (riga) immediatamente precedente che confmano con tale 
cluster elementare; poi in questo insieme S si identifica un cluster "vincitore" 
w che ha il maggior numero di pixel di confine con il cluster c e si uguaglia il 
suo numero identificativo con quello del cluster c. II risultato dell' operazione 
MERGEl sull'immagine binaria di Figura 15 e mostrato in Figura 16. 

I cluster deir insieme S diversi dal cluster vincitore vengono individuati 
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come cluster "perdenti" e ad essi e assegnata una rispettiva tripla 
comprendente il loro numero identificativo, il numero identificativo del 
rispettivo cluster vincitore e il numero di colonna cui appartengono (BOUND 
REPORTING, Report Processing). 

I numeri identificativi dei cluster perdenti vengono cambiati nei numeri 
identificativi dei rispettivi cluster vincitori mediante Toperazione MERGE2, il 
cui risultato e mostrato in Figura 17, A questo punto i cluster elementari 
appartenenti ad uno stesso oggetto sono individuati da uno stesso numero 
identificativo. 

Per comodita si possono assegnare ai diversi oggetti rappresentati 
numeri identificativi consecutivi con Poperazione di Cluster Sorting, 
come mostrato in Figura 18. 

Infme con Toperazione CLUSTER EXTRACTION, si estraggono 
dairimmagine originaria di Figura 19, associata aH'immagine binaria di 
Figura 15, i pixel corrispondenti a cluster elementari aventi lo stesso numero 
identificativo, isolando cosi i singoli oggetti rappresentati, come mostrato in 
Figura 20. 

Un possibile sistema di identificazione implementante il metodo 
descritto e mostrato in Figura 26. Esso comprende un sottosistema di 
scansione EXTRACTOR di pixel di un'immagine, un sottosistema di 
determinazione (FEATURES EXTRACTION, CLUSTERING CONDITION) di cluster 
elementari e un sottosistema di elaborazione CLUSTERING. II blocco circuitale 
Features Extraction ricava dai pixel scanditi i ripettivi tre antecedenti del 
criterio fiizzy (tono, 5, 5^), ® fomisce questa informazione al blocco in 
cascata Clustering Condition, che raggruppa i pixel in cluster elementari 
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applicando il criterio fiizzy. Infine il sottosistema di elaborazione 
Clustering dei cluster elementari produce in uscita pixels di oggetti 
rappresentati Cluster 1, CLUSTER N. 

II metodo di identificazione deirinvenzione pud essere 
convenientemente utilizzato per eseguire Toperazione di segmentazione 
(intra-spot segmentation) di spot di immagini rilevate da micro-array. Questa 
operazione consiste nel distinguere pixel di oggetti rappresentati in uno spot 
luminoso da pixel di sfondo o di rumore. 

Preferibilmente, ma non obbligatoriamente, airinterno dello spot sono 
identificate preliminarmente due zone: la zona True signal contenente pixel 
di segnale e la zona False signal contenente pixel di rumore o di sfondo. Per 
la separazione dei pixel dello spot in una zona contenente di pixel di segnale 
ed in una zona contenente di pixel di sfondo e utilizzato Talgoritmo fuzzy 
con i 3 antecedenti sopra citati, ciascun antecedente avente 3 funzioni 
membro (membership functions). 

Per I'apprendimento (learning) di un sistema fuzzy destinato ad 
implementare Talgoritmo, sono stati individuati gli spot piu rappresentativi 
della maggior parte dei casi che si possono verificare. Ad ogni pattern e 
associato un valore logico 1 o 0 a seconda che il pixel considerato sia un 
pixel di segnale o di sfondo. 

In Figura 27 e mostrato un diagramma di simulazione delFalgoritmo 
fuzzy in corrispondenza dei pattern di allenamento corrispondenti a due spot 
con livelli di grigio completamente difFerenti. 

Una volta completata I'operazione di segmentazione (intra-spot 
segmentation), ogni singolo spot del micro-array viene infine esaminato per 
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ottenere dei parametri caratteristici e indici di qualita. Questi parametri 
permettono di caratterizzare uno spot ed eventualmente di efFettuare anche 
un'analisi comparativa tra le immagini ottenute usando i coloranti Cy3 e Cy5 
per evidenziare Tibridazione in campioni di DNA. 

La seguente tabella illustra alcuni parametri utili per caratterizzare uno 
spot e il modo in cui essi sono calcolati, facendo riferimento alia figura citata. 



Parametro 



Centro dello spot 



Media del segnale 



Media del background 



Mediana del segnale 



Mediana del background 



FOMV 



Fattore di simmetria 



Area dello spot 



Area del background 



Area del segnale 



Fattore di correzione 



SigniFicato 



Coordinate del baricentro dello spot 
Figura 32 



Media dei livelli di grigio dei pixel dello spot 
appartenenti alia classe True signal 
Figura 33 



Media dei livelli di grigio dei pixel dello spot 

appartenenti al peri metro della shape 
Figura 34 



Mediana dei livelli di grigio dei pixel dello spot 

appartenenti classe True sigital 
Figura 3 5 



Mediana dei livelli di grigio dei pixel dello spot 
appartenenti al perimetro della shape 

Figura 36 



Media dell'uscita del sistema fiizzy per i pixel 
assegnati alia classe True signal. 



Rapporto tra altezza e larghezza del minimo 

rettangolo che include lo spot 
Figura 37 



Numero di pixel dello spot 
Figura 38 



Numero di pixel del perimetro della shape 
Figura 39 



Numero di pixel appartenenti alia classe True 
signal 
Figura 40 



DifFerenza tra la mediana del livello di grigio dei 
pixel del segnale e la mediana del livello di grigi^ ^ ^ 
dei pixel del background rii^^ 



Un parametro importante per la valutazione della qualita della selezionie<^ 
dei pixel estratti e quello indicato con la sigla FOMV, dato dalla media 
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deiruscita fiizzy per pixel di segnale. Quanto piu Tuscita fuzzy si approssima 
ad 1, tanto piii si puo considerare affidabile Tindividuazione di un pixel come 
pixel appartenente ad un oggetto rappresentato (pixel di segnale) e non come 
un pixel di sfondo. Facendo una media delPuscita fiizzy, si ottiene quindi un 
indice della qualita della selezione dei pixel di segnale tra i pixel di sfondo. 
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RIVENDICAZIONI 

1. Metodo di analisi di immagini rilevate da un micro-array sotto 
forma di uno o piu spot luminosi su uno sfondo, comprendente le seguenti 
operazioni: 

determinare forme e relative posizioni di detti spot sull'area 
sensibile di detto micro-array (array localization) generando una mappa 
binaria di pixels defmente confmi di detti spot luminosi su detto sfondo, 

isolare ciascuno di detti spot mediante un'operazione di estrazione 
(spot extraction) usando detta mappa binaria, 

esaminare ciascuno spot mediante un'operazione di segmentazione 
(intra-spot segmentation) identificando con un criterio prestabilito pixel 
appartenenti ad uno stesso cluster, 

per ciascuno di detti spot, determinare relativi parametri 
caratteristici ed indici di qualita, 

caratterizzato dal fatto che detta mappa binaria e generata con una 
tecnica di filtraggio morfologico comprendente le operazioni di: 

filtrare un'immagine rilevata da microarray con almeno un filtro 
morfologico generando una corrispondente immagine 'marker' dello sfondo 
(background); 

ricostruire detto sfondo (background) eflfettuando un'operazione di 
ricostruzione (reconstruction) su detta immagine 'marker', generando una 
corrispondente immagine ricostruita di sfondo; 

generare un'immagine filtrata della luminosita dello sfondo 
effettuando un'operazione top-hat su detta immagine ricostruita e detta 
immagine rilevata da microarray; 
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generare detta mappa binaria mediante un'operazione di 
thresholding su detta immagine fikrata della luminosita dello sfondo. 

2. II metodo della rivendicazione 1, in cui detta operazione di 
ricostruzione (reconstmction) e eseguita stabilendo preliminarmente che gli 
spot hanno forma circolare. 

3. II metodo della rivendicazione 1, comprendente inoltre 

filtrare il rumore corrompente detta mappa binaria mediante le 
seguenti operazioni: 

eseguire in successione due operazioni di erosion con maschere 
circolari di raggio diverso, 

eseguire un'operazione di dilation con una maschera circolare di 
diametro superiore al diametro massimo degli spot, generando una mappa 
binaria filtrata da rumore, 

utilizzare detta mappa binaria filtrata da rumore in detta 
operazione di estrazione (spot extraction). 

4. II metodo della rivendicazione 1, in cui gli spot di un' immagine 
'marker' sono generati con le seguenti operazioni: 

defmire un sistema di riferimento cartesiano per gli spot 
deir immagine rilevata da micro-array; 

eflfettuare in ordine le seguenti operazioni di filtraggio morfologico 
di detti spot con opening direzionali impieganti come elemento strutturante 
(structuring set) segmenti di lunghezza non superiore ad un massimo 
diametro degli spot, orientati rispettivamente 

lungo la bisettrice del primo e del terzo quadrante, 

lungo la bisettrice del secohdo e del quarto quadrante. 
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lungo Tasse delle ascisse, e 

lungo I'asse delle ordinate 
di detto sistema di riferimento cartesiano, gener£indo spot di una 
corrispondente immagine 'marker'. 

5. II metodo della rivendicazione 1, in cui detta operazione di 
estrazione (spot extraction) comprende le operazioni di: 

scandire (GENERAL CLUSTERING) pixel di un'immagine per 
colonna o per riga, associando a pixels adiacenti scanditi consecutivamente e 
corrispondenti a pixel della relativa mappa binaria aventi uno stesso valore 
logico attivo, una rispettiva quadrupla definente un cluster elementare 
costituita da un numero identificativo (INDEX), coordinate minima (xmin) e 
massima (xmax) e numero di colonna o di riga (y); 

individuare (MERGE 1) per ogni cluster elementare (INDEX=C) in 
una certa colonna o riga (i+1) un insieme di cluster elementari (S) nella 
colonna o riga immediatamente precedente (i) confmanti con detto cluster 
elementare (INDEX=C); 

individuare in detto insieme di cluster elementari (S) un cluster 
"vincitore" (INDEX=W) che ha il maggior numero di pixel di confine con 
detto cluster elementare (INDEX=C) e cluster "perdenti" diversi da detto 
cluster "vincitore", porre il numero identificativo di detto cluster elementare 
uguale al numero identificativo di detto cluster "vincitore" (W); 

porre (MERGE2) il numero identificativo di ciascuno dei clus4?>2r^^ 

Juj 

perdenti uguale al numero identificativo del rispettivo cluster vincitore (Wii'S^^ 
selezionare (CLUSTER SORTING, CLUSTER EXTRACTIOI^^^J:^ 



pixel rappresentanti spot luminosi estraendo dalF immagine originaria pixel di 
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cluster aventi lo stesso numero identificativo. 

6. II metodo della rivendicazione 1, in cui detto criterio prestabilito 
deH'operazione di segmentazione (intra-spot segmentation) consiste nei 
calcolare un valore caratteristico per pixel di uno spot con un algoritmo fiizzy 
contemplante le operazioni di: 

calcolare per detto spot il valore medio del tono di pixel di sfondo; 
detto algoritmo fiizzy avendo come antecedenti: 
il tono di un pixel, 

la distanza tra detto tono e il valore medio del tono dei rispettivi 
pixel di sfondo, 

il quadrato di detta distanza; 

giudicare detti pixel come appartenenti ad uno stesso cluster se 
detto valore caratteristico eccede una soglia prestabilita. 

7. II metodo della rivendicazione 6, comprendente le operazioni di 
definire con detto criterio prestabilito in ciascuno spot una prima zona (TRUE 
SIGNAL) contenente pixel di segnale ed una seconda zona (false signal) 
contenente pixel di sfondo o di rumore. 

8. II metodo della rivendicazione 6, in cui ciascun antecedente ha tre 
funzioni membro (membership). 

9. II metodo della rivendicazione 8, in cui dette funzioni membro 
(membership) sono funzioni gaussiane aventi media e varianza 
predeterminate. 

10. II metodo della rivendicazione 6, in cui detto algoritmo fiizzy ha 
cinque conseguenti. 

11. II metodo della rivendicazione 6, in cui detta operazione di 
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segmentazione (intra-spot segmentation) comprende: 

scandire (GENERAL CLUSTERING) pixel di un'immagine per 
colonna o per riga, associando a pixels adiacenti scanditi consecutivamente e 
soddisfacenti detto criterio prestabilito, una rispettiva quadrupla definente un 
cluster elementare costituita da un numero identificativo (INDEX), coordinate 
minima (xmin) e massima (x^ax) e numero di colonna o di riga (y); 

individuare (MERGEl) per ogni cluster elementare (lNDEX==C) in 
una certa colonna o riga (i+1) un insieme di cluster elementari (S) nella 
colonna o riga immediatamente precedente (i) confmanti con detto cluster 
elementare (INDEX=C); 

individuare in detto insieme di cluster elementari (S) un cluster 
"vincitore" (lNDEX=W) che ha il maggior numero di pixel di confine con 
detto cluster elementare (INDEX=C) e cluster "perdenti" diversi da detto 
cluster "vincitore", porre il numero identificativo di detto cluster elementare 
uguale al numero identificativo di detto cluster "vincitore" (W); 

porre (MERGE2) il numero identificativo di ciascuno dei cluster 
perdenti uguale al numero identificativo del rispettivo cluster vincitore (W); 

separare (CLUSTER SORTING, CLUSTER EXTRACTION) 
pixel di oggetti rappresentati in detto spot da pixel di sfondo estraendo pixel 
di cluster aventi lo stesso numero identificativo. 

12, II metodo della rivendicazione 7, comprendente calcolare per 
ciascuno spot parametri caratteristici ed indici di qualita scelti neirinsieme 
composto da: 

il valore medio del tono dei pixel di segnale di detta prima zona 
(TRUE SIGNAL); 
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le coordinate del baricentro dello spot; 

il valore medio del tono dei pixel perimetrici dello spot; 

il valore mediano dei toni dei pixel di segnale di detta prima zona 

(TRUE SIGNAL); 

il valore mediano dei toni di pixel perimetrici dello spot; 

il rapporto tra altezza e larghezza del minimo rettangolo 

circoscritto a detta prima zona (TRUE SIGNAL); 

il numero di pixel costituenti lo spot; 

il numero di pixel perimetrici dello spot; 

il numero di pixel di segnale di detta prima zona (TRUE SIGNAL); 
un fattore di normalizzazione del tono di pixel uguale alia 
differenza tra il valore mediano del tono dei pixel di detta prima 
zona (TRUE SIGNAL) e il valore mediano del tono dei pixel 
perimetrici dello spot; e 

il valore medio (FOMV) di detto valore caratteristico per i pixel di 
detta prima zona (TRUE SIGNAL). 



13. Dispositivo di analisi di immagini rilevate da un micro-array 
implementante il metodo della rivendicazione 1, comprendente un sistema di 
array localization configurato secondo un'architettura di rete neurale cellulare 
per il processamento di pixel di detto spot. 

14. II dispositivo della rivendicazione 13, avente un sistema di 
estrazione (spot extraction) di spot luminosi su uno sfondo di immagini 
rilevate da micro-array, comprendente 

un sottosistema di scansione (EXTRACTOR) di pixel di 
un'immagine; 
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un sottosistema di determinazione (Features Extraction, 
Clustering Condition) di cluster elementari formati da pixels adiacenti 
scanditi consecutivamente implementante Toperazione di estrazione (spot 
extraction) della rivendicazione 5, 

un sottosistema di elaborazione (CLUSTERING) di detti cluster 
elementari, producente in uscita pixels di cluster (CLUSTER 1, CLUSTER N) 
rappresentanti detti spot, 

15. II dispositivo della rivendicazione 14, avente un sistema di 
segmentazione (intra-spot segmentation) di spot luminosi su uno sfondo, 
comprendente 

un sottosistema di scansione (EXTRACTOR) di pixel di uno spot; 

un sottosistema fiizzy accoppiato al sottosistema di scansione, 
distinguente i pixel scanditi in pixel di segnale e pixel di sfondo o di rumore 
implementante il metodo della rivendicazione 6. 

16. II dispositivo della rivendicazione 15, in cui detto sottosistema 
fuzzy implementa Toperazione di segmentazione (intra-spot segmentation) 
della rivendicazione lie comprende 

un sottosistema di determinazione (FEATURES EXTRACTION, 
Clustering Condition) di cluster elementari formati da pixels adiacenti, 
costituente un sistema fuzzy, 

un sottosistema di elaborazione (CLUSTERING) di detti cluster / 

r 

elementari secondo il metodo della rivendicazione 1 1 , producente in uscita 
pixels di oggetti (CLUSTER 1, CLUSTER N) rappresentati in detto spot. 

17. Metodo di identificazione di pixel di un'immagine appartenenti ad 
uno stesso oggetto rappresentato su uno sfondo, comprendente 
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scandire pixel di detta immagine; 



calcolare un valore caratteristico per ciascun pixel con un 



aigoritmo fuzzy avente come antecedenti 

- il tono di un pixel, 

- la distanza tra detto tono e il valore medio del tono di rispettivi pixel di 
sfondo, 

- il quadrato di detta distanza; 
e contemplante le operazioni di: 



giudicare detti pixel come appartenenti ad uno stesso oggetto se detto valore 
caratteristico eccede una soglia prestabilita. 

18. II metodo della rivendicazione 17, in cui ciascun antecedente ha tre 
funzioni membro (membership). 

19. II metodo della rivendicazione 18, in cui dette funzioni membro 
(membership) sono funzioni gaussiane aventi media e varianza 
predeterminate. 

20. II metodo della rivendicazione 17, in cui detto aigoritmo fuzzy ha 
cinque conseguenti. 

21. Metodo di segmentazione (intra-spot segmentation) di spot 
luminosi su uno sfondo di un' immagine rilevata da un micro-array, 
comprendente Toperazione di esaminare ciascuno spot utilizzando il metodo 
della rivendicazione 17 per distinguere pixel di oggetti rappresentati in uno 
spot da pixel di sfondo o di rumore, 

22. II metodo della rivendicazione 17, comprendente: 



calcolare il valore medio del tono di pixel di sfondo; 



scandire (GENERAL CLUSTERING) pixel di un'immagine per 
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colonna o per riga, associando a pixels adiacenti scanditi consecutivamente e 
giudicati appartenenti ad uno stesso oggetto, una rispettiva quadmpla 
definente un cluster elementare costituita da un numero identificativo 
(index), coordinate minima (xmin) e massima (xmax) e numero di colonna o di 
riga (y); 



una certa colonna o riga (i+1) un insieme di cluster elementari (S) nella 
colonna o riga immediatamente precedente (i) confinanti con dettp cluster 
elementare (INDEX=C); 



"vincitore" (INDEX=W) che ha il maggior numero di pixel di confine con 
detto cluster elementare (INDEX=C) e cluster "perdenti" diversi da detto 
cluster "vincitore", porre il numero identificativo di detto cluster elementare 
uguale al numero identificativo di detto cluster "vincitore" (W); 

porre (MERGE2) il numero identificativo di ciascuno dei cluster 
perdenti uguale al numero identificativo del rispettivo cluster vincitore (W); 

separare (CLUSTER SORTING, CLUSTER EXTRACTION) 
pixel di oggetti rappresentati in detta immagine da pixel di sfondo estraendo 
pixel di cluster aventi lo stesso numero identificativo. 

23. II metodo della rivendicazione 21, in cui detta operazipne di 
segmentazione (intra-spot segmentation) comprende Toperazione di definire 
preliminarmente con detto criterio prestabilito in ciascuno spot una prima 
zona (true signal) contenente pixel di segnale ed una seconda zona (FALSE 
SIGNAL) contenente pixel di sfondo o di rumore. 

24. II metodo della rivendicazione 23, comprendente calcolare per 



individuare (MERGE 1) per ogni cluster elementare (INDEX=C) in 



individuare in detto insieme di cluster elementari (S) un cluster 
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ciascuno spot parametri caratteristici ed indici di qualita scelti neirinsieme 
composto da: 

il valore medio del tono dei pixel di segnale di detta prima zona 
(TRUE SIGNAL); 

le coordinate del baricentro dello spot; 

il valore medio del tono dei pixel perimetrici dello spot; 

il valore mediano dei toni dei pixel di segnale di detta prima zona 

(TRUE SIGNAL); 

il valore mediano dei toni di pixel perimetrici dello spot; 

il rapporto tra altezza e larghezza del minimo rettangolo 

circoscritto a detta prima zona (TRUE SIGNAL); 

il numero di pixel costituenti lo spot; 

il numero di pixel perimetrici dello spot; 

il numero di pixel di segnale di detta prima zona (TRUE SIGNAL); 
un fattore di correzione uguale alia differenza tra il valore mediano 
del tono dei pixel di detta prima zona (TRUE SIGNAL) e il valore 
mediano del tono dei pixel perimetrici dello spot; e 
il valore medio (FOMV) di detto valore caratteristico per i pixel di 
detta prima zona (TRUE signal). 
25. Metodo di analisi di immagini rilevate da un micro-array sotto 

forma di uno o piu spot luminosi su uno sfondo, comprendente le seguenti 

operazioni: 

determinare forme e relative posizioni di detti spot sull'area 
sensibile di detto micro-array (array localization) generando una mappa 
binaria di pixels definente confini di detti spot luminosi su detto sfondo. 
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isolare ciascuno di detti spot mediante un'operazione di estrazione 
(spot extraction) usando detta mappa binaria, 

esaminare ciascuno spot mediante un'operazione di segmentazione 
(intra-spot segmentation) utilizzando un criterio di similitudine prestabilito, 

per ciascuno di detti spot, determinare relativi parametri 
caratteristici ed indici di qualita, 

caratterizzato dal fatto che detta operazione di segmentazione (intra- 
spot segmentation) e eseguita secondo il metodo deila rivendicazione 21. 

26. Sistema di identificazione di pixel di un'immagine appartenenti ad 
uno stesso cluster rappresentato su uno sfondo, comprendente 

un sottosistema di scansione (EXTRACTOR) di pixel di immagine; 

un sottosistema fuzzy accoppiato al sottosistema di scansione ed 
identificante i pixel scanditi come pixel appartenenti ad uno stesso oggetto 
implementando il metodo della rivendicazione 17. 

27. II sistema della rivendicazione 26, in cui detto sottosistema fuzzy 
implementa il metodo della rivendicazione 22 e comprende 

un sottosistema di determinazione (FEATURES EXTRACTION, 
Clustering Condition) di cluster elementari formati da pixels adiacenti, 
costituente un sistema fuzzy, 

un sottosistema di elaborazione (CLUSTERING) di detti cluster 
elementari secondo il metodo della rivendicazione 22, producente in uscita 
pixels di oggetti (CLUSTER 1, CLUSTER N) rappresentati in d^t^ 
immagine. V^, 

28. II sistema deila rivendicazione 26, in cui detto sottosistema d 
determinazione (Features Extraction, Clustering Condition) di cluster 
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elementari costituisce un sistema fuzzy implementante il metodo 
segmentazione (intra-spot segmentation) della rivendicazione 21. 
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